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Circle separability queries in logarithmic time 

Greg Aloupis * Luis Barba * Stefan Langerman * 



Abstract 

Let P be a set of n points in the plane. In this paper we study a new variant of the 
circular separability problem in which a point set P is preprocessed so that one can quickly 
answer queries of the following form: Given a geometric object Q, report the minimum 
circle containing P and exluding Q. Our data structure can be constructed in 0(n log n) 
time using 0(n) space, and can be used to answer the query when Q is either a circle or a 
convex m-gon in O(logn) or 0(logn + logm) time, respectively. 
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^— i 1 Introduction 

The planar separability problem consists of constructing, if possible, a boundary that separates 
the plane into two components such that two given sets of geometric objects become isolated. 
Typically this boundary is a single curve such as a line, circle or simple polygon, meaning that 
each component of the plane is connected. Often there is an additional objective of minimizing 
some feature of the boundary, for instance the radius of a circle. 

Probably the most classic instance of this problem is to separate two given point sets with 



CZ) a circle (or a line, which is equivalent to an infinitely large circle). This problem was developed 

because of the several applications it has in pattern recognition and image processing [TJJ [Tj5] . 
Circle separability has been extensively studied and usually two variants are considered: The 
decision problem, where the existence of the separating circle is tested, and the Min-Max problem, 
where the minimum or maximum separating circle is to be determined. 

A separating line can be found, if it exists, using linear programing. In the plane this 
pq takes linear time by Megiddo's algorithm [T5]. For circle separability, O'Rourke, Kosaraju and 

Megiddo [TB] gave a linear-time algorithm for the decision problem (in fact spherical separability 
~f-! in any fixed dimension), improving earlier bounds [5J Q3]. They also gave an O(nlogn) algorithm 

for finding the largest separating circle and a linear time algorithm for finding the minimum 
separating circle between any two finite point sets. Extending these ideas, Boissonat et al. [6] 
gave a linear-time algorithm to report the smallest separating circle for two simple polygons, if 
£> it exists. 

The separability problem has also been studied when two point sets are to be isolated by 
the boundary of a simple polygon. Edelsbrunner and Preparata [TT] gave an 0(nlogn)-time 
algorithm to find a separating convex polygon with minimum number of vertices. This was 
shown to be optimal if the optimal separator has linear size. They also gave a O(kn) algorithm 
to find a convex polygonal separator with either k or k+1 vertices, where k is the size of the 
optimal solution. 

Aggarwal et al. PQ gave an 0(n log fc)-time algorithm to find the separating (convex) polygon 
with fewest vertices, between two nested convex polygons. Again, k is the size of the optimal 
separator. 
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Das and Joseph extended the problem to higher dimensions and proved that the problem of 
computing a separating polyhedron, having the minimum number of faces, for two nested convex 
polyhedra is NP-complete [7]. Suri and O'Rourke J20J gave a quadratic time algorithm to find a 
polygonal separator of two simple, not necessarily convex, nested polygons. This was improved to 
0(n log n) time by Wang [23] • Finally, Wang and Chan [23] gave an O(nlogn) algorithm to find 
a minimal polygonal separator between two arbitrary, not necessarily nested, simple polygons. 

In some situations, the given geometric objects and the separating boundary are constrained 
to lie within some subset of the plane, such as a simple polygon. For instance, Demaine et al. [8] 
studied the separation of point sets by chords and geodesic paths inside polygons. 

The online version of the separability problems, in which a preprocess is allowed to answer 
separability queries, is of special interest when a lot of queries are to be processed. This is the case 
in areas like geometric modeling involving objects in motion, where collision detection or largest 
empty space recognition queries are extensively used. Therefore, the online variants of some 
separability problems have been studied as well: Augustine et al [3J show how to preprocess a 
point set P, so that the largest circle, isolating P from a query point, can be found in logarithmic 
time. They also obtain this query time when P represents the boundary of a simple polygon. 

For the line separability problem, Edelsbrunner show that a point set P can be preprocessed 
in 0(n log n) time, so that a separating line between P and a query convex m-gon Q can be 
computed in 0(logn + logm) time [10]. In 3D space, Dobkin and Kirkpatrick show that two 
convex polyhedra can be preprocessed in linear time, so that a separating plane, if any exists, 
can be computed in 0(logn • logm) time, given any orientation and position of both polyhedra. 
In this case, n and m represent the size of each polyhedra. 

In this paper we show that a point set P on n points can be preprocessed in 0(n log n) time, 
using 0(n) space, so that for any given convex m-gon Q, we can find the smallest circle enclosing 
P and excluding Q in 0(log n + log m) time. This improves the 0(log n ■ log m) bound presented 
in [3], which is described in this paper as well. 

2 Preliminaries 

Let P be a set of n points in the plane and let Q be a convex m-gon. We say that every circle 
containing P is a P-circle. We also say that a circle C separates P from Q, or simply that C is 
a separating circle, if C is a P-circle and its interior does not intersect Q. Likewise, a separating 
line is a straight line leaving the interiors of P and Q in different halfplanes. 

Let C denote the minimum separating circle and let c' be its center. Note that C always 
passes through at least two points of P, since otherwise a smaller separating circle can always 
be found. In fact c' must lie on an edge of the farthest-point Voronoi diagram V(P), which is 
a tree with leaves at infinity [T7] . For each point p of P, let R(p) be the farthest-point Voronoi 
region associated with p. 

Let Cp be the minimum enclosing circle of P. If Cp is constrained by three points of P then 
its center, cp, is at a vertex of V(P). Otherwise Cp is constrained by exactly two points of P 
(forming its diameter), in which case cp is on the interior of an edge of V(P). If that is the case, 
we will insert cp into V(P) by splitting the edge where it belongs. Thus, we can think of V(P) 
as a rooted tree on cp. For any given point x on V(P) there is a unique path along V(P) joining 
Cp with x, throughout this paper we will denote this path by tt x . 

Given any point y in the plane, let C(y) be the minimum P-circle with center on y and let 
p{y) be the radius of C(y). Note that if y belongs to R(p) for some point p of P, then p{y) is 
given by d(y,p), where d(-, •) denotes the Euclidian distance between any two geometric objects. 
Finally, we say that y is a separating point if C(y) is a separating circle. 



3 Properties of the minimum separating circle 

In this section we describe some properties of C, and the relationship between c' and the 
farthest-point Voronoi diagram. These properties are not new. In fact most of the results in this 
section are either proved, stated, or assumed in [4j. 

Let CH(P) denote the convex hull of P. We assume that the interiors of Q and CH(P) are 
disjoint, otherwise there is no separating circle. Also, if Q and Cp have disjoint interiors, then 
Cp is trivially the minimum separating circle. 

The following useful property of the farthest-point Voronoi diagram was mentioned in [4j. 
We provide a short proof. 

Proposition 3.1. Let x be a point on V(P). The function p is monotonically increasing along 
the path n x starting at cp. 

Proof. Let u,v be two vertices of V(P). In [IS], they proved that if u is an ancestor of v, then 
p(u) < p(v). It only remains to prove that p is also monotonically increasing along every edge 
of n x . Let e be an edge of tt x contained in the bisector of the points p,p' of P, and consider a 
point y moving along e. Note that as d — d(y,p) — d(y,p') increases, the radius of C(y) also 
increases since y lies on R(p) fl R{p')- Thus, by moving the point y we obtain that p(y) increases 
monotonically on every edge along ir x . □ 

Observation 3.2. Every P-circle contained in a separating circle is also a separating circle. 

The following is stated in [4J . We provide a brief proof. 

Proposition 3.3. Let x and y be two points in M 2 . If z is a point contained in the segment 
[x,y], then C(z) C C{x) U C{y). 

Proof. Let a and b be the two points of intersection between C(x) and C(y). Let r = d(a, z) — 
d(b, z) and let C r {z) be the circle with center on z and radius r; see Figure [I] 

It is clear that C(x) n C(y) C C r (z) and since P C C(x) n C(y), we infer that C r (z) is a 
P-circle, therefore C(z) C C r (z). Furthermore, since C r (z) C C(x) U C(y) by construction, we 
conclude that C(z) C C{x) U C(y). D 



Observation |3 . 2| and Proposition 3.3 imply that if C(x) and C(y) are both separating circles, 
then for every z £ [x,y], C{z) is also a separating circle. This comes from the fact that C(z) 
is contained in C(x) U C(y) and Q lies outside of this union. Furthermore, this implies that 
the minimum separating circle is unique. Assume otherwise that C and C" are both minimum 
separating circles with centers c and c', respectively. Thus, for every point z in the open segment 
(c, c'), C(z) is also a separating circle contained in CU C . This means that C(z) has an smaller 
radius than C and C which would be a contradiction. 

The following is also demonstrated in [3], with a reworded proof. 

Lemma 3.4. Let x and y be two separating points on V(P). If z is the lowest common an- 
cestor of x and y in the rooted tree V(P), then C{z) is a separating circle; moreover p{z) < 
min{p(x),p(y)}. 



Proof. Since z is a common ancestor of x and y, we infer from Proposition 3.1 that p(z) < 
vo\n{p{x) 1 p(y)}. What remains is to prove that C(z) is a separating circle. Suppose that y ^ w x 
and x £ TT y , otherwise the result follows trivially since y or x will coincide with z. So, let e x 
{resp. e y ) be the edge incident to z, on the path joining z with x (resp. y). Consider the radial 
order of the edges incident to z lying between e x and e y . Let e' be the edge consecutive to e x in 




Figure 1: For every point z £ [x,y] the circle C(z) lies in the union of C(x) and C(y). 



this order. Therefore, e x ,e' and z all belong to the boundary of some Voronoi region R(p); see 
Figure [2j 

Let £ zp be the line through z and p. By the definition of V(P), £ ZyP intersects the boundary 
of R(p) only at the point z, and it separates x and y. Thus the intersection point z' between 



[x,y] and £ z ^ p belongs to R(p). Recall that Proposition 3.3 implies that C(z') is a separating 



circle. Trivially, C{z) C C(z') since the latter is obtained by expanding C(z) while anchoring 



it to p. Since C{z) is a P-circle contained inside a separating circle, by Observation 3.2 C(z) is 
also a separating circle. □ 




C(z') 



Figure 2: Illustration for Lemma pO] 



The following is also proved in [4]. 



Theorem 3.5. Let s be a point on V(P). If s is a separating point, then c' belongs to tt s . 
Proof. Proceed by contradiction and assume that c' does not belong to ir s . Let z be the lowest 



common ancestor of c' and s. Note that z lies in tt s and therefore z/c'. In this case, Lemma 3.4 



in conjunction with Proposition 3.1 imply that C{z) is a separating circle with p(z) < p(c') which 
is a contradiction. □ 

Given a separating point s, we claim that if we move a point y continuously from s towards 
cp on n s , then C(y) will shrink and approach Q, becoming tangent to it for the first time when 
y reaches c'. To prove this claim in Lemma |3.8[ we introduce the following notation. 

Let re be a point lying on an edge e of V(P) such that e lies on the bisector of p,p' € P. Let 
C~(x) and C + (x) be the two closed convex regions obtained by splitting the disk C(x) with the 
segment [p,p'}. Assume that x is contained in C~(x). 

Observation 3.6. Let x,y be two points lying on an edge e ofV(P). If p{x) > p{y), then 
C+{x) C C+{y) and C-(y) C C~(x). 




Figure 3: The result of Observation 3.6 when p(x) > p(y) 



Using this simple observation that can be seen in Figure [3j we obtain the following general- 
ization. 

Proposition 3.7. Let s be a point on V(P) and let x and y be two points on n s . If p(x) > p{y), 
then C+{x) C C+{y) and C~{y) C C~(x). 



Proof. Note that if x and y lie on the same edge, then the result holds by Observation 3.6 If they 
are on different edges, we consider the path $ = (x, vq, . . . , Vk, y) contained in tt s joining x and y, 



such that Vi is a vertex of V(P), i £ {0, . . . , fc}. Thus, Observation 3.6 and Proposition 3.1 imply 



that C+{x) C C+{v ) C . . . C C+(v k ) C C+{y) and that C~(y) C C~{v k ) C . . . C C-(« ) C 
C-(x). D 

Note that C = C(c') must intersect the boundary of Q. Otherwise, we could always push 
c' closer to the root on V(P), while keeping it as a separating point. Furthermore, since Q is 



convex and C contains no point of Q in its interior, the intersection consists on only one point. 
From now on we refer to <f>' as the tangency point between C and Q. 

We claim that 4>' lies on the boundary of C + (c'). Assume to the opposite that <f>' lies on 
C~(c'). Let e > and let c £ be the point obtained by moving c' a e distance towards to cp 



on V(P). Note that by Proposition 3.1 p(c e ) < p(c'). In addition, Proposition 3.7 implies that 
C~(c e ) C C~(c'). Since we assumed that <f>' lies on the boundary of C~{d), we conclude that </>' 
does not belong to C(c £ ). This implies that, for e sufficiently small, C(c e ) is a separating circle 
which is a contradiction to the minimality of C. 

Lemma 3.8. Let s be a separating point. If x is a point lying on it s , then C{x) is a separating 
circle if and only if p(x) > p(c'). Moreover, C is the only separating circle that intersects Q. 

Proof. Let s be a separating point. We know by Theorem |3 . 5 1 that c' belongs to tt s . Let x\ and 



x 2 be two points on 7r s such that p{x\) < p(c') and p(c') < p(x 2 ). Proposition 3.7 implies that 
C + (c') C C + (xi) and since 4>' belongs to the boundary of C + (c'), we conclude C(x\) contains 
<f>' in its interior. Therefore C{x\) is not a separating circle. 

On the other hand, C(x 2 ) contains no point of Q. Otherwise, let q €E Q be a point lying in 
C(x 2 ). Two cases arise: Either q belongs to C~(x 2 ) or q belongs to C + (x2). In the former case, 
since p{s) > p(x 2 ), q G C~(x 2 ) C C~(s) — a contradiction since C(s) is a separating circle. In 



the latter case, since p(x 2 ) > p(c'), Proposition 3.7 would imply that q belongs to the interior of 



C which would also be a contradiction. □ 

4 The algorithm 

The basis of our algorithm is to find a separating point s and from there, perform a binary search 
on tt s to find a separating circle tangent to Q with center on this path. This is a rather intuitive 
solution, which was also used in 0]. However we use some additional geometric properties to 
reduce the time complexity. 

4.1 Preprocessing 

We first compute V(P) and cp in O(nlogn) time [TS]. V(P) can be stored as a binary tree with 
n (unbounded) leaves, so that every edge and every vertex of the tree has a set of pointers to the 
vertices of P defining it. Every Voronoi region is stored as a convex polygon and every vertex p 
of P has a pointer to R(p). If cp is not a vertex of V(P), we split the edge that it belongs to. 

We want our data structure to support binary search queries on any possible path tt s of V(P). 
Thus, to guide the binary search we would like to have an oracle that answers queries of the 
following form: Given a vertex v of 7r s , decide if c' lies either between cp and v or between v 



and s in n s . By Lemma 3.8 we only need to decide if C(y) is a separating circle. 

We will use an operation on the vertices of V(P) called FindPointBetween with the fol- 
lowing properties. Given two vertices u, v in 7r s , FindPointBetween(u, v) returns a vertex 
z that splits the path on 7r s joining u and v into two subpaths. Moreover, if we use our or- 
acle to discard one of the subpaths and to proceed recursively on the other, then we want 
FindPointBetween(u, v) to guarantee that this recursive process ends after O(logn) steps. 
That is, after O(logn) iterations, the search interval becomes only an edge of ir s containing c'. 

A data structure build on top of V(P) that support this operation was presented in |18j . 
This data structure can be constructed in 0(n) time and uses linear space by storing a constant 
number of pointers on each vertex of V(P). 



4.2 Searching for c' on the tree 

Recall that if Cp is a separating circle then it is a trivial solution. Since Q is a convex m-gon, this 
can be checked easily in O (log to) time |10) . Thus we will assume that Cp is not the minimum 
separating circle, which implies that Cp intersects Q. 

To determine the position of c' on V(P), we first find a separating point s and then search 
for c' on 7r s using our data structure. 

To find s, we construct a separating line L between P and Q. This can be done in 0(logn + 
log to) time [TO]. Let p L be the point of P closest to L and assume that no other point in P 
lies at the same distance; otherwise rotate L slightly. Let L± be the perpendicular to L that 
contains p L and let s be the intersection of Lj_ with the boundary of R(p L ) ; see Figure [IJ We 
know that Lj_ intersects R{p L ) because L can be considered as a P-circle, containing only p L , 
with center at infinity on L±. 




Figure 4: Construction of s. 



Since s is on the boundary of R(p L ), C(s) passes through p L . Furthermore C(s) is contained 
in the same halfplane defined by L that contains P. So C(s) is a separating circle. Assume that 
s lies on the edge xy of V(P) with p(x) > p(y) and let 7r s = (uq — s, u\ — y, . . . , u r — cp) be the 
path of length r + 1 joining s with cp in V(P). Theorem 3.5 implies that c' lies on ir s . 



It is possible to use our data structure to perform a binary search on the vertices of tt s , com- 
puting, at each vertex v, the distance to Q and the radius of C(v). This way we could determine 
if C(v) is a separating (or intersecting) circle. However, this approach involves computing the 
distance to Q at each step in O(logm) time, and thus takes (9(logn • logm) time. This was the 
algorithm given in 0]. 

It is worth noting that if the distance to our query object can be computed in O(l) time, then 
the described algorithm takes O(logn) time. This is the case when instead of being an m-gon, 
Q is cither a circle or a point. The next result follows. 

Observation 4.1. After preprocessing a set P of n points in O(nlogn) time, the minimum 
separating circle between P and any given circle or point can be found in O(logn) time. 

In the case of Q being a convex TO-gon, an improvement from the 0(lognlogTO.) time algo- 
rithm can be obtained by strongly using the convexity of Q. 



To determine if some point v on ir s is a separating point, it is not always necessary to compute 
the distance between v and Q. One can first test, in constant time, if C(v) intersects a separating 
line tangent to Q. If it does not intersect it, then C(v) would be a separating circle and we can 
proceed with the binary search. Otherwise, we can try to compute a new separating line tangent 
to Q not intersecting C(v). The advantage of this is that while doing so, we reduce the portion 
of Q that we need to consider in the future. This is done as follows. 

Compute the two internal tangents L, V between the convex hull of P and Q in 0(logn + 
log m) time. The techniques to construct these tangents are shown in Chapter 4 of [21] . 

Let q and q' be the respective tangency points of L and L' with the boundary of Q. We will 
consider the clockwise polygonal chain ip= [q = q , . . . ,q/, — q'] joining q and q' as in Figure pj 




Figure 5: The construction of ip. 



Recall that <j>' denotes the intersection point between C and the boundary of Q and note 
that the tangent line to C at <$>' is a separating line. Therefore, <j)' must lie on an edge of ip since 
no separating line passes through any other boundary point of Q. 

If q = q' , then <f>' = q and thus, we can forget about Q and compute the minimum separating 



circle between P and q. As mentioned previously, by Observation 4.1 this takes O(logn) time. 

Assume from now on that q ^ q' , as shown in Figure pj For each edge e, = qiqi+i (0 < 
i < k — 1) of tp, let li be the line extending that edge. By construction, we know that each £i 
separates P and Q. We say that a point x on li but not in e, lies to the left of e, if it is closer 
to qt, or to the right if it is closer to <fo+i. 

Our algorithm will essentially perform two parallel binary searches, the first one on ir s and 
the second one on (p, such that at each step we discard either a section of 7r s or a linear fraction 
of if. As we search on it s , every time we find a separating circle, we move towards cp. When we 
confirm that a P-circle intersects Q, we move away from cp. 

As mentioned, we attempt to confirm if the current vertex v being analyzed corresponds to 
a separating point. To do that, we compare C(v) to some separating line £i for intersection, i.e., 
in constant time. If C{v) is a separating circle, we instantly discard the section of the path lying 
below v on V(P). If C(v) does intersect £i, we make a quick attempt to check if C(v) intersects 
Q by comparing C(v) and the edge e^ for intersection. If so, v is not a separating point and we 
can proceed with the binary search on tt s . Otherwise, the intersection of C(v) with £i lies either 
to the left or to the right of e,. However, in this case we are not able to quickly conclude whether 



C(v) intersects Q or not. Thus, we temporarily suspend the binary search on V(P) and focus 
on Civ), using it to eliminate half of ip. Specifically, the fact that C(v) intersects ti to one side 
of e, (right or left) tells us that no future P-circle on our search will intersect ti to the other 
side of e, . This implicitly discards half of ip from future consideration, and is discussed in more 
detail in the Theorem that follows. Thus, in constant time, we manage to remove a section of 
the path n s , or half of ip, which will give the desired time bound. The entire process is detailed 
in Algorithm [l] 

Algorithm 1 Given 7r s = (uq = S,U\ = y, . . . ,u r — cp) and ip = [q = qa, . . . , qi~ = q'], find the 
edge of 7r s containing c' 
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Define the endpoints of the subpath of n s containing c', u <— s, v <— cp 
Define the initial search interval on ip, a <— 0, b <— k 
if u and v are consecutive vertices and b = a + 1 then 

End and return the segment S — [u,v] and the segment H = [q a ,qb] 
end if 

Let z •(— FindPointBetween(m, v), j <— L^J 
Let ej ^— qjqj^i and let £j be the line extending ej 
if b > a + 1 then 

Compute p(z) and let S <— d(z, £j), A ^— d(z, ej) 
else 

Compute p(z) and let 5 4— d(z, ej), A 4— d{z, ej) 
end if 

if p(z) < S, that is C(z) is a separating circle then 
Move forward on 7r s , u <— z and return to step [3] 
else 

if p(z) > A , that is if C(z) is not a separating circle then 

Move backward on tt s , v <— z and return to step [3] 
else 

if C(z) intersects lj to the left of ej then 

We discard the polygonal chain to the right of ej, b 4— max{j, a + 1} 
else 

We discard the polygonal chain to the left of ej, a -s— j 
end if 

Return to step [3] 
end if 
end if 



Theorem 4.2. Algorithmu] finds the edge of tt s containing c' in 0(logn + logm) time. 

Proof. Our algorithm maintains two invariants. The first is that C(u) is never a separating circle 
and C(v) is always a separating circle. To begin with, C(u) = C(s) is a separating circle, while 
C(v) = Cp is not a separating circle. If either of these assumptions does not hold, the problem is 



solved trivially, without resorting to this algorithm. Changes to u and v occur in steps 14 or 17 
and in both cases the invariant is preserved. Thus, c' always lies on the path joining u with v. 

The second invariant is that <$>' , the tangency point between C and Q, always lies on the 
clockwise path joining q a with qb along ip. We already explained that the invariant holds when 
a = and b = k, corresponding to the inner tangents supporting P and Q. Thus we only need 



to look at steps 20 and[22{ where a and b are redefined. 



We will analyze step 20 however step 22 is analogous. In step 20 we know that C(z) intersects 
£j to the left of ej and that ej does not intersect C(z). We claim that for every point w lying on 
an edge of ir s , if C(w) is a separating circle that intersects £j, then it intersect it to the left of 
ej. Note that if our claim is true, we can forget about the polygonal chain lying to the right of 
ej since no separating circle will intersect it. To prove our claim, suppose that there is a point w 
on 7r s , such that C(w) is a separating circle and C(w) intersects tj to the right of ej. Let x and 
x' be two points on the intersection of lj with C(w) and C(z), respectively. Suppose first that 



p(w) < p(z) and recall that by Proposition 3.7 since x' lies on C + (z) C C + (w), x' lies in C(w). 
Thus, both x and x' belong to C(w) which by convexity implies that ej is contained in C(w). 
Therefore C(w) is not a separating circle which is a contradiction. Analogously, if p(w) > p(z), 
then ej is contained in C{z) which is directly a contradiction since we assumed the opposite; our 
claim holds. 

Note that in each iteration of the algorithm, a, b, u or v are redefined so that either a linear 
fraction of ip is discarded, or a part of n s is discarded and a new call to FindPointBetween is 
performed. Recall that our data structure guarantees that O(logn) calls to FindPointBetween 
are sufficient to reduce the search interval in n s to an edge [TS] . Thus, the algorithm finishes in 
0(log n + log to) iterations. 

One extra detail needs to be considered when b = a+l. In this case only one edge e = [q a , q a +i] 
remains from <p, and <j>' lies on that edge. Thus, if the line £ extending e intersects C(z) but 
e does not, then either step [20] or [22] is executed. However, nothing will change in these steps 
and the algorithm will loop. In order to avoid that, we check in step [5] if only one edge e of (p 
remains. If this is the case, we know by our invariant that <p' belongs to that edge and therefore 
we continue the search computing the distance to e instead of computing the distance to the 
line extending it. This way, the search on ip stops but it continues on ir s until the edge of V(P) 
containing c' is found. 

Since we ensured that every edge in V(P) has pointers to the points in P that defined it, 
every step in the algorithm can be executed in 0{\) time. Thus, we conclude that Algorithm IT] 
finishes in 0(logn + log to) time. 

Since both invariants are preserved during the execution, Lemma |3.8| implies that the algo- 
rithm returns segments [u,v] from ir s containing c', and [<7 a ,<7b] from ip containing <fi'. □ 

From the output of Algorithm [TJ it is trivial to obtain c' in constant time, so we conclude the 
following. 

Corollary 4.3. After preprocessing a set P of n points in Oinlogn) time, the minimum sepa- 
rating circle between P and any query convex m-gon can be found in 0(log n + log to) time. 
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